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ABSTRACT 

We present and analyse integral-field observations of six type-II QSOs with z=0.3-0.4, 
selected from the Sloan Digital Sky Survey (SDSS). Two of our sample are found to be 
surrounded by a nebula of warm ionized gas, with the largest nebula extending across 
8" (40 kpc). Some regions of the extended nebulae show kinematics that are consistent 
with gravitational motion, while other regions show relatively perturbed kinematics: 
velocity shifts and line widths too large to be readily explained by gravitational mo- 
tion. We propose that a ~ 20 kpc x20 kpc outflow is present in one of the galaxies. 
Possible mechanisms for triggering the outflow are discussed. In this object, we also 
find evidence for ionization both by shocks and the radiation field of the AGN. 
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1 INTRODUCTION 

Powerful active galactic nuclei present unique opportuni- 
ties for investigating the formation and evolution of mas- 
sive galaxies, as well as the interplay between the active 
galactic nucleus (AGN) and its environment. Type-II (i.e., 
obscured) radio-quiet quasars (QSOs hereinafter) are of par- 
ticular interest for this. Their high luminosities should afford 
their detection across large ranges in redshift, while their 
lack of powerful radio jets should simplify the interpreta- 
tion of the observable properties of their host galaxies and 
environments. Moreover, most of the growth of super mas- 
sive black holes occurs in obscured quasars (e.g. Martinez- 
Sansigre et al. 2005). In addition, unlike their Type-I (un- 
obscured) counterparts, the optical emission from the host 
galaxy is unlikely to be overwhelmed by the AGN. 

The existence of Type-II QSOs was required by unifica- 
tion models of active galaxies (e.g. Barthel 1989; Antonucci 
1993 and references therein) and, although long sought-after, 
it was only with the Sloan Digital Sky Survery (SDSS) that 
these objects were detected in significant numbers (Zakam- 
ska et al. 2003). To date, there have been a number of spec- 
troscopic investigations of Type-II quasars (e.g. Zakamska 
et al. 2003; Villar-Martm et al. 2008), but these have been 



relatively superficial and based predominatly on data from 
the SDSS, focussing primarily on classifying the objects, un- 
derstanding the nature of the powering mechanism and the 
obscuring structure, and testing unification models. 

While these studies have provided a wealth of impor- 
tant information, little work has been carried out to char- 
acterize the gaseous properties of Type-II quasars. In this 
letter we present integral field spectroscopic observations 
for 6 moderate redshift type II QSOs, 5 of which are radio 
quiet (Liaghz < 10 31 ergs" 1 Hz" 1 Sr" 1 ). Our main goals 
are to investigate the possible existence of extended ionized 
nebulae around the galaxies, and characterize for the first 
time their kinematic, morphological and excitation proper- 
ties. Throughout this letter, we assume a flat universe with 
#o=71 km s" 1 Mpc~\ fi A =0.73 and fi m =0.27. 



2 OBSERVATIONS AND REDUCTION 

The sample we use in this investigation consists of one radio- 
loud and five radio-quiet type-II quasars (QSOs) at z~0.3- 
0.4. This redshift range was used to enable us to observe 
from [NeV]A3426 to [OIII] A5007, and also to minimise any 
redshift evolution effects. All of our sources were selected 
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from the sample of type-II QSOs originally identified by Za- 
kamska et al. (2003) in the SDSS. The objects were selected 
to have relatively high [OIII] A5007 emission line fluxes, i.e, 
) 1 x 10 -15 erg s _1 cm" 2 , in order to maximise the proba- 
bility of detecting extended nebulae. We also required that 
source coordinates permit observations from Calar Alto. 

Our observations were taken on 2008 January 8 and 
10, and on 2008 July 1 and 2, at the 3.5m telescope of the 
Calar Alto observatory, Spain, using the Potsdam Multi- 
Aperture Spectrograph (PMAS: Roth et al. 2005) in its lens- 
array mode. The atmospheric conditions were variable. Dur- 
ing the January observing run the humidity was high (i.e., 
non-photometric). In the July run the conditions were pho- 
tometric. The seeing was relatively stable, ranging between 
1.2" and 1.6" (or ~6-8 kpc). The PMAS lens-array com- 
prises a bundle of 256 fibres, coupled to a square array of 
256 lenses, and resulted in contiguous sampling of the sky. 
The optics were adjusted to give a spatial scale of l"/spaxel 
(lens) and a field of view of 16" x 16" . Cross-talk between ad- 
jacent fibres is estimated in less than a 3% when using a pure 
aperture extraction (see below). We used the V300 grating, 
covering a spectral range of ~4250-7650A, wherein a num- 
ber of interesting emission lines lie. The full-width at half 
maximum (FWHM) of the instrumental profile, measured 
from the strong sky line [OI] A5577A, was 7.4A ±0.5A. Our 
observations are summarized in Table 1. 

Data reduction was performed using R3D (Sanchez 
2006), in combination with IRAF packages and E3D 
(Sanchez 2004). The reduction consists of the standard steps 
for fibre-based integral-field spectroscopy (Sanchez 2006), 
including bias and flat-field corrections, tracing and extrac- 
tion of spectra, geometrical distortion correction, and cali- 
bration of the wavelength and flux scales. Finally, the spec- 
tra were assembled into a data-cube. We made no correc- 
tion for galactic extinction, because it is low for all sources; 
uncertainties introduced by not applying the correction are 
expected to be substantially smaller than those introduced 
by the flux calibration. 



3 ANALYSES 

3.1 Emission line images 

In order to examine the spatial distribution of the line emit- 
ting gas associated with the galaxies, we have reconstructed 
images of the [OIII] emission doublet and, where the signal- 
to- noise ratio was sufficient, [Oil]. First, we summed the 
line flux along the dispersion axis of the data cube, to pro- 
duce line plus continuum images. For each combination of 
object and emission line, we used a fixed wavelength range 
selected to include the wings of the line, and also to al- 
low for any velocity shifts between lenses, corresponding to 
an image bandwidth of between llOA and 220A. In the 
case of the [OIII] doublet, we included both the 4959A line 
and the 5007A line. Next, we reconstructed images for the 
continuum emission from a spectral region nearby to the 
emission line, and with the same wavelength bandwidth as 
the relevant line plus continuum image. These were sub- 
tracted from the line plus continuum images, to produce 
pure line images. The reconstructed emission line images are 
shown in Fig. 1. In each image, the lowest contour level rep- 
resents the 3cr detection limit. We do not detect spatially 
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Figure 1. Reconstructed emission line images for our sample. 
North is at the top and East is to the left. The axes give the spatial 
offset, in arcsec, from the continuum peak. The contours are plot- 
ted for 1, 2, 4, 6.7, 13, 25 ,50 X 3<r when possible. The line emission 
from SDSS J084041. 08+383819 and SDSS J155059.30+395029.5 
is spatially extended on scales of several tens of kpc; their signif- 
icant spatial variations in kinematic properties and in line ratios 
confirms their spatial extensions (Figs. 2 and 5). In the rest of 
the sample, we find no evidence spatially extended line emission 
at the resolution of our observations (~6-8 kpc). 



extended line emission from SDSS J092014. 11+453157.3, 
SDSS J162545.11+310632.2, SDSS J172603.07+602115.6 or 
SDSS J173938.64+544208.6 (see column 8 of Table 1). This 
is predominantly due to the faintness (or lack of) extended 
line emission from these galaxies. We find no clear trend 
between detection/non-detection of extended nebulae and 
the properties of the emission form the central spaxel or 
the radio properties (Table 1; c.f. Christensen et al. 2006; 
Husemann et al. 2008). The four galaxies without detec- 
tions of extended nebulae will not be considered further 
in this letter. For SDSS J155059. 30+395029.5 and SDSS 
J084041. 08+383819, we also show [Oil] /[OIII] maps in Figs. 
2 and 5, respectively. 



3.2 Velocity and FWHM maps 

Maps of the FWHM and line of sight velocity were pro- 
duced for [OIII] A5007, typically the brightest emission line 
in these spectra. We fitted the spectral region around [OIII] 
A5007 with a single Gaussian velocity profile atop a linear 
continuum, making use of MPFIT and other routines from 
the Markwardt IDL library. We resticted ourselves to fit 
only those lenses where the [OIII] A5007 velocity profile had 
a peak signal-to-noise ratio of ^7. Although the kinemat- 
ics are complex in some lenses, we have nevertheless fitted 
single Gaussians because our goal is obtain a general char- 
acterisation of the kinematic properties. FWHM values have 
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Table 1. Summary of the observations. Column 3 lists the radio type: RQ=radio quiet (l>\AGHz < 10 31 ergs~ 1 H z -1 Sr~ 1 ); 
RL=radio loud (L\AGHz ~& 10 31 ergs -1 H z~ l Sr~ 1 ). Column 7 gives the la sensitivity of the observations relative to that of SDSS 
J084041. 08+383819. 8. Column 8 gives the mean surface brightness of the [OIII] A5007 line, in units of 10 — 16 erg s — 1 cm -2 arcsec - 2 , at a 
radius of 3"from the peak of the optical continuum. Columns 9 lists the FWHM of[OIII] A5007 after deconvolution from the instrumental 
profile (units of km s ), in the spaxel corresponding to the continuum flux maximum. 



Target 


z 


Type 


Date 


Exp. (s) 


Notes 


a/a 




FWHM 


SDSS J084041.08+383819.8 


0.313 


RQ 


10/01/08 


8x1800 


high humidity 


1.0 


1.6±0.1 


1330±130 


SDSS J092014.11+453157.3 


0.402 


RQ 


08/01/08 


2x1800 


cloudy 


2.0 


^0.6 


700±100 








10/01/08 


2x1800 


cloudy 








SDSS J155059.30+395029.5 


0.347 


HQ 


01/07/08 


6x1200 


photometric 


1.8 


4.8±0.2 


920±90 


SDSS J162545.11+310632.2 


0.379 


RL 


02/07/08 


6x1200 


clear 


3.6 


^0.9 


540±40 


SDSS J172603.07+602115.6 


0.333 


RQ 


01/07/08 


6x1200 


photometric 


2.2 


^0.3 


900±100 


SDSS J173938.64+544208.6 


0.384 


HQ 


02/07/08 


7x1200 


clear 


2.8 


^0.2 


660±70 



been corrected for instrumental broadening. FWHM and Av 
have typical la uncertainties of ~100 km s~ . 



4 RESULTS AND DISCUSSION 
4.1 SDSS J155059. 30+395029. 5 

Although this source meets our criteria for being radio-quiet, 
it does emit a significant radio flux density. The flux den- 
sity measurement from the FIRST survey (Becker, White 
& Helfand 1995), which has a spatial resolution of 5" (24 
kpc), is 10.6+0.5 mjy. From the NVSS survey (Condon et 
al. 1998), with a resolution of 45" (220 kpc), the flux density 
is measured to be 13.1+0.6 mjy. The fact that the NVSS 
measurement is significantly higher suggests that there is 
emission on large spatial scales, and this is almost certainly 
emission from radio jets. 

This object is associated with a spatially extended emis- 
sion line nebula (Fig 1), with an integrated [OIII] A5007 lu- 
minosity of 1.5±0.2xl0 43 erg s" 1 . The [OIII] AA4959,5007 
and [Oil] A3727 doublets extend across 8" (40 kpc) and 6" (30 
kpc), respectively. H/3 and [NeV] A3426 are also spatially 
extended. The nebula is elongated approximately along the 
North-South axis, and shows an excess of line emission on 
East side of the emission peak relative to the West. 

The kinematic properties vary significantly across the 
nebula: 130^FWHM<1070 km s _1 and -20«SAv<310 km 
s _1 (Fig. 2). The outer regions of the nebula have narrower 
emission lines, i.e. ,^420 km s - , consistent with gravita- 
tional motion, and the blue wing is weak or undetected 
therefrom. We also note that regions with lower FWHM 
tend to have higher line of sight velocities (Fig. 4). 

The highest FWHM are measured in the spaxel at the 
spatial zero, and in several spaxels to the immediate North 
and East thereof (up to 1070 km s _1 ). In these spaxels, a 
substantial blue wing is also apparent in the velocity profiles 
(Fig. 4), with a total detected extent of 4"x4"(20 kpc x 20 
kpc). Since these kinematic properties are not centred on 
the flux peak of the line and continuum, they are unlikely 
to be the result of contamination by the central regions of 
the AGN. We define a 'detection' of the blue wing as fol- 
lows. We perform two separate fits to the line profile: one 
using a single Gaussian and another using two Gaussians. 
We formally 'detect' the blue wing when both of the follow- 
ing conditions are met: (1) the two-Gaussian fit shows an 



excess of flux on the blue side of the line profile, and (2) the 
reduced \ 2 ls ^16 higher for the two-Gaussian fit than it is 
for the single Gaussian fit. 

The presence of these highly perturbed kinematics sug- 
gests the existence of an outflow (e.g. Jarvis et al. 2003). In 
the absence of a detection of photospheric absorption lines, 
we will assume that the quiescent gas is at the systemic ve- 
locity, and that the blue wing is due to an outflow. Next, 
we attempt to quantify the observed properties of this out- 
flow. After summing the spaxels from which we detected 
the blue wing, we fitted the line profile of [OIII] A5007 using 
two Gaussians (Fig. 4). We use [OIII] A5007 because it is 
the highest signal to noise line in this spectrum. From the 
resulting fit, we take the relatively broader, blueshifted com- 
ponent to represent the outflowing gas: it has a (convolved) 
FWHM=1007 km s -1 , is blueshifted from our fiducial veloc- 
ity zero by u =285 km s _1 , and has a [OIII] A5007 flux of 
1.7 xl0~ 14 erg s _1 cm~ 2 . Scaling this to R/3 using the [OIII] 
A5007 / H/3 ratio measured in the same aperture, we obtain 
a H/3 flux of 2.7 x 10 -15 erg s _1 cm -2 , which corresponds to 
Lh/j=1.1 x 10 42 erg s _1 . While this parameterisation is likely 
to be a gross oversimplification of reality, it is a necessary 
one given our limited information about the true geometry 
and dynamics of the outflow. Assuming spherical outflow, 
the gas outflow rate is then 

3v L H prn v 



M = 



ff 



mhuBpctHp 



(1) 



where m p is the proton mass, r is the radius of the 
outflow, n is the density of the outflowing gas, hvufi is the 
energy of an H/3 photon and is the effective H/3 re- 

combination coefficient (e.g. Dopita & Sutherland 2003). 
For the outflow radius, we adopt the maximum projected 
spatial extent of the blue wing from the optical nucleus 
(^2"or ^3.0 xlO 22 cm). In addition, we assume n^lOOO 
cm~ 3 , based on the strong detection of the forbidden [Oil] 
doublet; the [SII] AA6716,6731 doublet profile also suggests 
densities lower than 1000 cm -3 . Hence, we obtain a very 
approximate outflow rate of > 2M©yr -1 . 

A galactic scale ouflow, as observed in relatively nearby 
galaxies, is one possible physical explanation. Such outflows 
can have spatial scales as large as ~10 kpc, outflow rates of 
~ 1-100 Moyr -1 , and are powered by active star formation 
with rates of >10 M0yr _1 (e.g. Heckman, Armus & Miley 
1990; Lehnert & Heckman 1996; Heckman et al. 2000). If 
the entire FIRST flux density were due to star-formation, it 
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would imply a star formation rate in excess of 1000 Mqyt , 
which would be more than sufficient to power the outflow 
from SDSS J155059.30+395029.5. 

In addition, we note that kinematically perturbed emis- 
sion lines have been detected from the extended ionized 
gas of radio-loud quasars and of powerful radio galax- 
ies, and are generally attributed to the radio jet ac- 
tivity (e.g. Clark et al. 1998). If the radio emission 
from SDSS J155059.30+395029.5 is dominated by emis- 
sion from the AGN, rather than powered by star forma- 
tion, could such a radio source power the outflow in SDSS 
J155059. 30+395029. 5? Assuming a radio spectral index of 1, 
and using the NVSS radio flux density, we calculate a radio 
luminosity of ~4xl0 31 erg s _1 Hz _1 Sr" 1 . Then, using the 
relation L ki „ = W 28±2 L° r ^± - 05 from Birzan et al. (2004), 
we estimate the kinetic jet power to be between 2xl0 37 and 
lxlO 45 erg s _1 . For an outflow velocity v a = 285 km s~\ 
a gas density of 1 cm' 3 and injection time of 10 Myr, the 
required energy injection rate is ~ 10 43 erg s" 1 (see e.g. 
equation 2 in Nesvadba et al. 2006), in agreement with the 
higher end of the range in possible kinetic jet power: thus, 
the radio source does appear able to power the ionized gas 
outflow. 

A further possibility is that the outflow is powered by 
the radiation field of the active nucleus: an outlfow is pre- 
dicted for quasars accreting near to, or above, the Eddington 
limit (e.g. King & Pounds 2003). In principle, outflows of this 
kind can escape from the central region of the AGN (King 
& Pounds 2003; Pounds et al. 2003a,b) and their spatial ex- 
tent can be of the order of tens of kpc (Silk & Rees 1998; 
King 2003). Assuming the relationship between the bolo- 
metric and [OIII] A5007 luminosities L& ; ~ 1000 x Ltoni] 
(e.g. Zakamska et al. 2003; Richards et al. 2006), the ob- 
served [OIII] luminosity of L[oiii] =1.5±0.2x 10 43 erg s -1 
implies Lj, ; ~ 1.5 x 10 46 erg s _1 , which in turn implies 
a black hole mass of 1.2 xlO 8 M0. Assuming a typical ra- 
diative efficiency of 5 per cent (Martinez-Sansigre & Taylor 
2009), this implies an accretion rate of ~ 5M0yr _1 . King & 
Pounds (2003) argue that the outflow rate will be roughly 
the same as the accretion rate, in which case the outflow rate 
would be in agreement with the rate of > 2M0yr~ we es- 
timated for SDSS J155059.30+395029.5. We note that this 
kind of outflow might provide a plausible explanation for 
the perturbed kinematics in some of the extended nebulae 
associated with radio galaxies and quasars, which heretofore 
have been assumed to be caused by their powerful radio jets. 

The [Oil] /[OIII] ratio varies substantially across the 
nebula: 0.3<[OII]/[OIII]<1.7. At x=l"N (5 kpc N), y=2"W 
(10 kpc W), this ratio reaches its maximum, suggesting a 
relatively low state of ionization. In addition, the [OIII]/H/3 
ratio of 2.1±0.8 is also consistent with low a ionization state 
at this position (Fig 3). The detection of spatially extended 
[NeV] A3426 suggests that the hard radiation field of the ac- 
tive nucleus contributes significantly to the ionization of the 
nebula. Interestingly, [Nelll] A3869 / [OIII] A5007=l. 0+0.3, 
which is unusually high and is consistent with ionization by 
shocks, but is ten times higher than would be expected for 
photoionized gas (e.g. Dopita & Sutherland 1996). We sug- 
gest that shock-ionization is a consequence of the spatially- 
extended outflow in this source. 

Moreover, in the spaxels where the lines are broadest, 
we find that the FWHM of H/3 is systematically higher than 
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Figure 3. Spectrum from the lens at x=l" (5 kpc) N, y=2" 
(10 kpc) W in SDSS J155059.30+395029.5. Notice the relative 
weakness of [OIII] A5007 compared to [Oil], [Nelll] and H/3; 
[OIII]A5007/H/3=2. 1+0.8 and [NeIII]A3869/[OIII]A5007=1.0±0.3 
suggest ionization by shocks. 

those of [OIII] AA4959,5007 by between 100±40 and 220±40 
km s _1 , suggesting that the blue wing / outflow has a lower 
ionization state than the kinematically quiescent gas. This is 
also consistent with the effects of shocks (Clark et al. 1998). 

4.2 SDSS J084041. 08+383819. 8 

We measure spatial extents of 5" (23 kpc) and 6" (27 kpc) 
for the [OIII] AA4959,5007 and [Oil] A3727 emission lines, 
respectively. The total luminosity of the [OIII] A5007 line 
is 9.8+1.5 xlO 41 erg s" 1 . While the [Oil] flux peaks at the 
position of the optical continuum peak, the [OIII] emission 
flux reaches its maximum l"(5 kpc) W,l"(5 kpc) S of this 
position. The spatial distribution of the line emission is elon- 
gated towards the West. 

Fig. 5 shows the FWHM, line of sight velocity, and [Oil] 
A3727 / [OIII] AA4959,5007 ratio maps. FWHM ranges be- 
tween 580 and 1150 km s _1 , with the highest values be- 
ing measured at, or adjacent to, the spatial zero. A promi- 
nent red wing is visible in the [OIII] velocity profile in the 
3"x3"(14 kpc x 14 kpc) centred on the spatial zero; it is 
not clear whether these kinematics originate in the extended 
gas, or are instead due to contamination by emission from 
the central regions of the AGN. Line of sight velocity varies 
between -80 and 230 km s _1 . The relatively broad emission 
lines at some spatial positions reveal the presence of highly 
perturbed gas in the nuclear region, or close to it. The [Oil] 
A3727 / [OIII] AA4959,5007 ratio varies from 0.5-1.1 across 
the nebula, with lower values closer to the spatial peak of 
the [OIII] emission. Overall, the line ratios are consistent 
with ionization by shocks or the radiation field of an AGN. 



5 SUMMARY AND DISCUSSION 

In this letter we have presented new integral-field ob- 
servations of 6 type-II QSOs at 0.3^z^0.4, obtained 
with PMAS on the 3.5m telescope at Calar Alto. Two 
of these sources, SDSS J155059.30+395029.5 and SDSS 
J084041. 08+383819.8, were found to be associated with ex- 
tended ionized nebulae, which have minimum spatial extents 
of 27 kpc and 40 kpc, respectively. Both nebulae exhibit var- 
ied kinematic properties: some regions show relatively quies- 
cent kinematics, consistent with gravitational motion, while 
in other regions of the nebulae the kinematics are highly 
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Figure 2. Results from the datacube of SDSS J155059.30+395029.5. From left to right: the [Oil] A3727 / [OIII] AA4959,5007 map; the 
[OIII] A5007 FWHM map (km s _1 ); the map of the velocity shift of [OIII] A5007 (km s _1 ), where we adopt the velocity at the spatial 
zero as our fiducial velocity zero. The range in [Oil] A3727 / [OIII] AA4959,5007,FWHM and velocity shift is typical of the spatially 
extended ionized gas associated with powerful active galaxies (see text). The total intensity map is similar to that shown in Fig. 2. 
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Figure 5. Results from the datacube of SDSS J084041. 08+383819.8. From left to right: the [Oil] A3727 / [OIII] AA4959,5007 map; the 
[OIII] A5007 FWHM map (km s _1 ); map of the velocity shift of [OIII] A5007 (km s —1 ), adopting the velocity at the spatial zero as our 
fiducial velocity zero. The range in [Oil] A3727 / [OIII] AA4959,5007,FWHM and velocity shift is typical of the spatially extended ionized 
gas associated with powerful active galaxies (see text). The total intensity map is similar to that shown in Fig. 1. 



perturbed, with FWHM exceeding 1000 km s 1 . We inter- 
pret the kinematic pattern in SDSS J155059.30+395029.5 
as an outflow of ~285 km s _1 velocity which extends across 
an area of 20 kpc x 20 kpc, with a mass ejection rate of 
> 2Af0yr _1 . Plausible mechanisms to power this apparent 
outflow include a starburst, the radio jets, or the AGN ra- 
diation field. 

Furthermore, the nebula associated with SDSS 
J155059.30+395029.5 shows evidence for AGN photoioniza- 
tion, based on the detection of extended [NeV] emission, and 
evidence for ionization by shocks, based on a relatively high 
[NeIII]/[OIII] ratio (see e.g. Dopita & Sutherland 1996). 

Overall, the luminosities, sizes, excitation properties an 
kinematics of the two extended nebulae are consistent with 
those of radio galaxies and type I QSOs at similar redshifts 
(see e.g. Baum et al. 1988; Willott et al. 1999). We suggest 
that AGN radiation powered outflows or star formation, ei- 
ther of which are able to explain the outflow from SDSS 
J155059. 30+395029. 5, might also be able to explain the sim- 
ilar kinematic properties shown by some nebulae associated 
with radio galaxies and quasars - properties which are often 
assumed to be powered by radio jets. 
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Figure 4. Left: Velocity profiles of the [OIII] AA 4959,5007 dou- 
blet in 1550+39, to illustrate the trend between FWHM and line 
of sight velocity Av. In this figure, an arbitrary flux scale was 
used. The 'broad' aperture is the sum of the 6 fibres showing the 
highest FWHM, i.e., [-2:0,0:1]. The 'narrow' aperture corresponds 
to 6 fibres in the Northern region of the nebula, i.e., [-1:0,2:3], and 
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'intermediate' is the sum of fibres [-1:0,-1]. The values for Av and 
FWHM are given in km s . Right: An illustration of how we 
'detect' blue wings. Here we show the observed [OIII] A5007 line 
in the 'broad' aperture (solid smooth line) fitted using two Gaus- 
sians (dashed red lines). The reduced \ 2 f° r this fit is 18.7, which 
is substantially lower than the value 94.2 obtained from a single 
Gaussian fit (not shown). In both panels, the flux density scale is 
in arbitrary units. 
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